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Abstract. A U (3) model proposed by Iachello for superconductiv- 
ity in cuprate materials is analyzed. The model consists of s and d 
pairs (approximated as bosons) in a two-dimensional Fermi system 
with a surface. The transition occurs between a phase in which the 
system is a condensate of one of the bosons, and a phase which is a 
mixture of two types of bosons. In the current work we have inves- 
tigated the validity of the Bogoliubov approximation, and we used a 
reduced Hamiltonian to determine a phase diagram, the symmetry of 
the phases and the temperature dependence of the heat capacity. 

For more than a decade there has been a debate regarding the na- 
ture of the symmetry of the macroscopic wavefunction for the copper 
oxide superconductors. According to Miiller bulk sensitive experi- 
ments support substantial s symmetry, while surface sensitive exper- 
iments yield d symmetry for the macroscopic wave function. Miiller 
proposed [1J that a reconciliation of the conflicting experiments was 
possible if the superconducting wavefunction were a sum of two com- 
ponents, namely, one of s symmetry and one of d symmetry, and var- 
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ied as a function of the distance from the surface. 

Following that work by Miiller, a theoretical framework for 
cuprate superconductors was proposed by Iachello [2J, based on the 
analogy with atomic nuclei. In the present work we report our analysis 
of phase transitions based on Iachello's model. Also we have inves- 
tigated temperature dependence of the heat capacity. We propose a 
macroscopic wavefunction for realizing space dependence. 

In Eq{T]we give Iachello's Hamiltonian for 2N spin j fermions 
on a plane. This is a boson Hamiltonian Q , which consists of three 
types of bosons s, d + and <i_. 

H = e s U i A+e d (£j + + JLJj\+uo(tff*S£\+ (1) 

/ aJ. y>J. A\ A\ /Vj- /S4- A\ A\ \ j I A\!~ a\!~ a\ A, \ 

+ u 2 ( d+d + d+d+ + d_d_d-d- J + u 2 ( d + d_d+d- j + 
+ vq ^d^_d+ +dLd-j y^f) + v 2 [d^dLss + $ s 1 d+d-^j 
Introducing the operators Q 



■ ^j. a. ^j. a. 

h s = s's hd+ = d\_d+ hd- = d_d- 

fid = h d+ + h d _ I = h d+ -h d - N = h d+ + h d - 

the Hamiltonian can be rewritten as 



(2) 



H = £ s n s + e d n d + u n s (n s - 1) + (3) 
+ 2 U 2"d (nd - 2) + +^u 2 n 2 d + -(u 2 - u 2 ) F + 
+ von d n s + v 2 (^d^cLLss + $ $ d+d-^J 
Its algebraic structure is U (3) 

i AJL. A. A-L. A. A, A. ^J- A 

h s = s's hd = d\_d + + d_d- I = d+d + — d_d— 



Q + = d l + d- f_ = -dls+s*d + f + = d\s-$d- 
Q = dld+ R+ = d\s + s^d- R = dls + s Jf d+ 

This algebra contains a subalgebra 50(3) formed from the genera- 
tors T + ,T , I , and a Casimir operator defined as (note a factor of 2 
difference from Ref. Q) 

A. '"J A\ A\ A\ A\ /V) 

w 2 = r+Ti + Tiri+z 2 
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Following Iachello this Hamiltonian can be split into an SO (3) 
symmetric part, which will contain I and W 2 terms, and a U (1) Cg> 

A A 

U (2) symmetric part generated by n s and l 1 n dl Q± 1 respectively. 
When the Hamiltonian has SO (3) symmetry the system is in phase 77, 
which corresponds to a mixture of s- and d- bosons. When the Hamil- 
tonian has U (1) (g> U (2) symmetry the system is in phase I, which 
corresponds to either an s-boson condensate (phase la) or a J-boson 
condensate (phase Ib)\2\. 

As shown in Ref. [2J an analytical solution could be obtained 
either for phase I or phase II, but in intermediate situations one could 
only obtain an approximate or a numerical solution. Since we have 
three kinds of bosons in the model and a macroscopic number of parti- 
cles, we may use the Bogoliubov approximation [4J. This approxima- 
tion consists of replacing in the Hamiltonian (HJ) the operators which 
create and annihilate bosons in the condensed state by a c-number and 
neglecting all terms higher than degree two in operators which anni- 
hilate and create bosons in other states. After this replacement the 
Hamiltonian becomes the bilinear form 

H d = H + (£ s + voN d -n)n s + (4) 

H = (£ d -^)N d + ^(2u 2 + u 2 )Nj + ^(2u2-u 2 )l 2 

where d±,d± — > y^N±, N d = N++N-, I = N + - N_ and jl is a chem- 
ical potential. The chemical potential is determined as usual from the 
condition on the mean total number of particles. N d is determined 
by variation. To diagonalize the bilinear form @ we introduce the 
squeezing transformation 

[/, = exp|(^_^. 

We compare the Bogoliubov approximation with the exact numerical 
solution in Fig [T] The right column of Fig [l] shows the occupancy 
of the ^-bosons. We observe different regions of the parameter cor- 
responding to different phases of the system. The left column shows 
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Figure 1. Relative error in the ground state energy (left y-axis, dashed line), 
calculated using the Bogoliubov approximation, and exact mean s-bosons 
population (right y-axis, solid line) versus parameters of the Hamiltonian. 



The fixed parameters are N = 100, / = 0, e s 
v = 1.5, v 2 = 1.3 



Ed = 10, Uq = U2 



0, 



that the Bogoliubov approximation gives remarkably good results ei- 
ther when the system is in phase la or phase II. The failure of the 
Bogoliubov approximation in phase lb is expected since the s-level is 
not macroscopically occupied in that phase. (Another Bogoliubov ap- 
proximation based on macroscopic occupancy of only J-bosons could 
be constructed for phase lb) 

In order to investigate the transition between phase I and phase 
77 the original Hamiltonian © can be reduced to the form 



H = £n d -j-W 2 



(5) 



by setting £ s = vq = 2u2 = —2u' 2 = — V2, uq = and e = £ d + ^V2- 
This Hamiltonian contains the essential features required to describe 
the transition between the phases. The system is in phase I for V2 = 
and it is in phase II for e = 0. In the first case the system reduces to a 
simple two level system of N non-interacting bosons, and in the sec- 
ond case we have pure pairing interaction. In order to have a quantita- 
tive measure of s-wave and d-wave symmetry we define two fractional 
weight operators 



. 1 . 
gd = j^n d 



1 . 
N 
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Figure 2. S-wave symmetry measure of the ground state calculated for vari- 
ous values of £ and V2- Lighter tones denote higher values. 

and the fractional weight measures of a state l^ 1 ) 

Sy=(V\§ s \V) Dy = (*¥\§ d \V) 
with the properties 

D = 1, D N/2 = ^, D N = 
1 

So = 0, S N / 2 = -, Sn = 1 

where the subscripts 0, f and N indicate |¥) = \0,N) , | f , f ) , \N, 0) 
respectively. 

For pure s-wave symmetry S^> and yield one and zero, re- 
spectively, and vice- versa for pure d- wave symmetry. For mixed sym- 
metries, S and D vary between zero and one. 

A phase diagram of the system described by the reduced Hamil- 
tonian © may be obtained readily by calculating the s-wave or the 
J-wave symmetry measure in the exact ground state for different val- 
ues of V2 and e. 

As seen from Fig.|2]the system has three distinct phases. When 
the system is in phase / the number of s- and d- bosons is definite; 
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when it is in phase II the number of s- and J-bosons is indefinite. The 
transition between the phases can be controlled by changing £ or vi- 

The strength of the scattering interaction (V2) could be con- 
trolled by the doping of the superconducting sample. It is known 
that doping affects electron-electron and electron-phonon scattering 
0, therefore the desired scattering on the Fermi surface could be 
achieved by varying the doping of the sample (Meir Weger, private 
communication). 

Having the exact numerical solution we can examine the ther- 
modynamical properties of this model by calculating the free energy 
and the density matrix of the system at inverse temperature /3 

Z = Tre~ pH F = -4lnZ p = -e~ pH 

j8 H Z 

We can then explore the behavior of the derivatives of the free energy 
with respect to the temperature. The first and the second derivatives of 
the free energy are shown to be continuous functions at the explored 
domain of parameters. However, there is a jump in the second deriva- 
tive when the temperature scale is of the same order as the difference 
between the ground and first excited states (Fig [21). 

This jump does not become a singularity point when the number 
of particles tends to infinity, since the width of the peak does not scale 
as N~ a , and therefore there is no second order phase transition at this 
point|6|. 

Although there is no second order phase transition, changing the 
temperature does produce a transition between phases having different 
symmetries. If the system is in phase I at T = 0, we can calculate the 
mean ensemble s-wave and rf-wave symmetry measures as a function 
of the temperature 

S = Tr(s,p) D = Tv{g d p). 

In Fig. |4] we present the temperature dependence of S and D. 
We observe that both symmetry measures approach j when the tem- 
perature is increased, which corresponds to phase 77 of the system. 
If the system is initially in phase II, it would remain in this phase 



7 




100 200 300 400 500 600 



Temperature 

Figure 3. The entropy (solid line) 
and the specific heat of the system 
(dashed line) as a function of the 
temperature. For T = the system 
is in phase lb. 




N=100 
e = v 2 = - 1 
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Figure 4. Mean symmetry mea- 
sures as functions of the temper- 
ature. The solid line denotes the 
s-wave symmetry measure and the 
dashed line denotes the d-wave 
symmetry measure. For T = the 
system is in phase lb. 



for any temperature. It is tempting to propose the temperature as the 
control parameter for the transition, but some things should be kept 
in mind. First, the temperature scale of the transition may be higher 
than the critical temperature of the superconductor. This can only be 
known after the model is fit to experimental data and the values of the 
Hamiltonian parameters © are obtained. Second, if we would like the 
temperature to be the control parameter of the transition, there must 
be a temperature gradient within the copper-oxide plane of the su- 
perconductor, so that the bulk of the superconductor will have higher 
temperature than its surface. This is not likely to happen since the su- 
perconductor is in thermal equilibrium and the typical surface depth 
is of the order of only several atomic layers. 

To sum up, in this work we have used a theoretical model pro- 
posed by Iachello Q for a transition from d-wave to mixed d- and 
s-wave symmetry. A phase diagram of the model was obtained, which 
showed a distinct separation between the three phases. Doping of the 
superconducting sample and temperature were proposed as possible 
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control parameters for the phases of the system. 
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